Protracted Screening in the Periodic Anderson Model 
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The asymmetric infinite-dimensional periodic Anderson model is examined with a quantum Monte 
Carlo simulation. For small conduction band filling, we find a severe reduction in the Kondo scale, 
compared to the impurity value, as well as protracted spin screening consistent with some recent 
controversial photoemission experiments. The Kondo screening drives a ferromagnetic transition 
when the conduction band is quarter-filled and both the RKKY and superexchange favor antiferro- 
magnetism. We also find RKKY-driven ferromagnetic and antiferromagnetic transitions. 
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Introduction Kondo lattice materials, stoichiometric 
systems generally with U or Ce atoms in the valence 
shell, have been studied intensely over the last few 
decades. They display a wide variety of behaviors in- 
cluding the well-known transport anomalies associated 
with a strongly enhanced electronic mass, as well as mag- 
netic (antiferromagnetic, incommensurate and ferromag- 
netic), paramagnetic, and superconducting ground states 
The magnetic nature of the ground state is deter- 
mined by the strength of the hybridization between the 
f-electrons with the delocalized band states ^ . If the hy- 
bridization is small, then the RKKY exchange dominates 
and the ground state is magnetic. If the hybridization 
is large, then either the Kondo screening removes the f- 
moments, or charge fluctuations delocalize the f-electrons 
destroying their moments, and the ground state is a Pauli 
paramagnet. 

The periodic Anderson model (PAM) is thought to de- 
scribe the magnetic and transport properties of these 
materials. Both it and the single impurity Anderson 
model (SIAM) have been extensively studied; neverthe- 
less, the detailed nature of the orbitally non-degenerate 
PAM phase diagram as well as the differences in the 
Kondo screening between the two models are unknown. 
In this Letter, we present the first rigorous calculations of 
the phase diagram of the asymmetric PAM in infinite di- 
mensions. We find antiferromagnetic, ferromagnetic, and 
paramagnetic ground states, depending upon conduction 
band filling rid and model parameters. In the paramag- 
netic state, we find that the temperature dependence of 
the Kondo screening retains the same qualitative features 
as the SIAM, whereas the quantitative features are quite 
different. Specifically, when Tid « [rid w n/ « 1) the 
Kondo scale is strongly suppressed (enhanced) and the 
temperature dependence of the screening is protracted 
(contracted) compared to the SIAM. These differences 
between the PAM and SIAM cannot be removed by 
rescaling and may provide insight into recent controver- 
sial photoemission experiments 0J. 

The PAM Hamiltonian is, 
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In (1), d(f)^J destroys (creates) a d(f)-electron with spin 
a on site i on a hypercubic lattice, U is the screened 
Coulomb repulsion potential for the localized f-states 
and V characterizes the mixing between the two sub- 
systems. This model retains the complications of the im- 
purity problem, which is known to be non-perturbative, 
including moment formation and screening and is further 
complicated by interactions between the moments due to 
RKKY and super-exchange mechanisms. 

Formalism A simplifying method which allows for an 
exact solution of the lattice problem in a non-trivial limit 
is necessary. Such a method was proposed by Metzner 
and Vollhardt Q| who observed that the irreducible self- 
energy and vertex-functions become purely local as the 
coordination number of the lattice increases. As a con- 
sequence, the solution of most interacting lattice mod- 
els may be mapped onto the solution of a local corre- 
lated impurity coupled to an effective bath that is self- 
consistently determined [pi. We refer the reader to the 
above references and recent reviews for further details on 
the method J(|. 

To solve the remaining effective impurity problem, we 
use the Anderson impurity algorithm of Hirsch and Fye 
0. In order to model the Ce-based Kondo lattice ma- 
terials, we place the correlated f-band below the Fermi 
level (so rif « 1) and adjust the conduction band filling 
by varying the Fermi level. Thus, beginning at [3 = 10, 
we choose e/ and so that rif — 1 and rid is the de- 
sired value. When the temperature is changed, we keep 
et — ed fixed and vary et to conserve the total number 
of electrons. For the results presented here the variation 
of nf from one was less than a few percent, and statisti- 
cal error bars are less than 7%. Note that our results are 
symmetric under rid —> 2 — rid, and the PAM becomes 
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more correlated for small values of rid and 2 — n<j. 

Results Fig. 1 shows the Kondo scale for the PAM 
and the SIAM versus d-band filling when the f-band is 
half-filled. The Kondo scales are obtained by extrapola- 
tion Xim P (T — > 0) = I /To, where x imp (T) is the addi- 
tional local susceptibility due to the introduction of the 
effective impurity into a host of d-electrons We see 
that at the symmetric limit (nj — nd = 1) the Kondo 
scale for the PAM is enhanced compared to the Kondo 
scale of SIAM as has been found earlier However, 
far from the symmetric limit the Kondo scale for the 
PAM is strongly suppressed. 

The main consequence of this suppression is that the 
temperature dependence of the screening in the PAM is 
quite different from that of the SIAM. This is shown 
in Fig. 2 where the screened local f- moments Txff for 
both models are plotted versus temperature. As expected 
from the impurity problem Txff displays a log- linear 
T-dependence in some temperature interval. However, 
concomitant with the differences in To, the screening of 
the PAM and the SIAM local moments are quite differ- 
ent: for most values of rid these two curves cannot be 
made to overlap by rescaling their temperature dependen- 
cies (nor is it possible, in the PAM case, to make the 
curves of different d-fillings overlap by rescaling T). At 
high temperatures, the screened local moments are iden- 
tical for the PAM and the SIAM, and as the temperature 
is lowered below Tqsiam the screening begins; however, if 
rid <SC ri/ i n d ~ nf = 1), the screening is significantly pro- 
tracted (contracted) in temperature for the PAM com- 
pared to that in the SIAM. (The contracted screening 
region, as well as the enhancement of the Kondo scale, 
when nd = rif = 1, depends strongly upon the correlation 
energy U, and diminishes when U is small; i.e., for the 
less correlated data set in Figs. 1 and 2). This behavior 
is consistent with recent controversial photoemission ex- 
periments on single-crystals of Ce-based heavy-fermion 
compounds, where the temperature-dependence of the 
screening peak is much slower than what would be ex- 
pected for a SIAM f|. 

The Drude weight, D, calculated by extrapolation of 
the current-current correlation function po| , is shown in 
the inset to Fig. 1. For all of the data shown in Fig. 1, D 
is quite small and the effective electron mass (not shown) 
is large m* /m — l/Z > 15, where Z is the quasiparticle 
renormalization factor. In the symmetric limit, where 
a gap opens in the single-particle density of states, we 
have D = Z = 0. However, consistent with what is seen 
in the Kondo scale for the PAM, D and Z also become 
very small when When the d- filling falls, there 

are no longer enough d-electrons to completely screen 
the f-moments. Thus, Tq falls quickly as the d-band is 
doped away from half filling. If this process continues, 
the d-band becomes depopulated, and then only acts as 
a hopping path between the f-levels. In the limit, nd — > 0, 
the PAM may be mapped onto a strongly correlated sym- 



metric Hubbard model pdf | (with an on-site correlation 
U and a strongly reduced hybridization) , which is known 
to open a gap in the single-particle density of states and 
have Z = and D = |(|. Thus, when rid = or n<j = 1, 
we find that both the f and the d-band densities of states 
vanish at the Fermi surface. However, for all other values 
of nd explored, the f and d-density of states remain finite 
(with the f-DOS only moderately enhanced) indicating 
that the system remains metallic. 

It is clear from the Drude weight shown in the inset 
to Fig. 1, that the unscreened moments have a dramatic 
effect on the Fermi-liquid properties of the system. This 
may also be seen by examining the electronic distribu- 
tion function n(e k ) = T^ n [G dd (e fc , iui n ) + G ff (e k , iui n )] 
where G dd and are the fully dressed d- and f-band 
Green's functions calculated with the QMC. dn(ek) / dek, 
is calculated by numerically evaluating the derivative of 
the above sum. The width of this distribution at low 
temperatures (shown in Fig. 3) gives an estimate of the 
single-particle scattering rate. It must go to zero if a 
Fermi-liquid is to form. This appears to happen when 
rid = 0.8,0.6 and 0.4; however, for rid — 0.2, it is not 
clear whether a Fermi-liquid forms. When nd ^$ 0.4, 
there is a protracted region in T of strong spin-flip scat- 
tering, beginning at T > TqsiaMi an d extending down 
to very low temperatures. 

In fact, due to magnetic ordering for nd ^$ 0.6, we 
find no compelling evidence for non Fermi liquid ground 
states. Fig. 4. shows the magnetic phase-diagram of 
the half-filled f-band PAM. T c is obtained by extrap- 
olating or interpolating the magnetic susceptibility as- 
suming the mean-field form \ ^ (t-t ) • There are 
two well-known exchange mechanisms that are usually 
responsible for the formation of these magnetic ground- 
states. The superexchange which results from the ex- 
change of the local f-electrons via the hybridization with 
the d-band; this exchange always favors antiferromag- 
netic order of the half-filled f-band and becomes strongly 
suppressed as the d-band is filled towards rid — 1 p2| - 
The Ruderman-Kittel-Kasuya-Yosida (RKKY) exchange 
which results from the consequent scattering of a d- 
electron from two f-momcnts; this exchange varies (in 
sign and magnitude) as a function of the d-band fill- 
ing (ljjl^]. The inset to Fig. 4 shows the difference 
between the staggered and uniform RKKY exchange, 
AJrkky = Jrkky{<\ = Q) - J rkky (q = 0), where 
Q is the wave- vector corresponding to the corner of the 
first Brillouin zone and 

JRKKYW = -7^ J fd E G,M ( k > i"n)G dd (k + q, iu n ) . 
n,k 

(2) 

Here, Jfd — —8V 2 /U, is the effective exchange between 
the f and d-band H|. AJrkky would be proportional 
to the energy gained by the formation of a ferromagnet 
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versus an antiferromagnet in the f-band if RKKY were 
the only exchange mechanism present. 

In the symmetric limit, n/ = = 1, for the more 
strongly correlated model (larger U), there is an insta- 
bility towards antiferromagnetism as expected from the 
sign of the RKKY-exchange. (The superexchange, which 
always favors antiferromagnetism, is expected to be neg- 
ligible in this limit |l2|). For the less correlated model 
(smaller U) the Kondo scale is larger (Fig. 1) and the 
d-electrons effectively screen away the local moments be- 
fore a magnetic transition can occur H, so there is no 
antiferromagnetism in these cases. 

As the d-filling decreases from 1, the transition remains 
commensurate and T c drops quickly to zero. However, 
as the system becomes more correlated, the size of this 
antifcrromagnetic region increases. As decreases fur- 
ther, the system remains a paramagnetic Fermi liquid 
until the d-band approaches quarter-filling. Here, both 
the superexchange and the RKKY-exchange still favor 
antiferromagnetism; nevertheless, the system has a fer- 
romagnetic transition. The rather low transition temper- 
ature (compared to T c at the symmetric limit), which is 
of the order of the Kondo scale, indicates that the mecha- 
nism behind this ferromagnetism might be related to the 
Kondo screening effect. 

To see how Kondo screening can lead to a ferromag- 
netic exchange, assume the size of the Kondo polarization 
cloud is 1/kp p| so the screening cloud is almost lo- 
cal. In this situation, every time a d-electron hops to the 
neighboring d-level it breaks its resonance with the local 
f-electron and regains the lost resonance-energy only if 
the f-electron on the neighboring site has the same spin 
orientation and the d-level is unoccupied. 

This Kondo-exchange mechanism is substantiated by 
examining the staggered charge susceptibility of the d- 
levels (not shown) which is strongly enhanced at low tem- 
peratures near quarter-filling. Hence, each site occupied 
by a d-electron tends to have its neighboring sites unoc- 
cupied, optimizing the above mentioned mechanism. In 
addition, if the quarter-filled d-electrons actually form a 
spin-polarized staggered charge-density wave, where the 
spins of the d-electrons are aligned opposite to the ferro- 
magnetically ordered f-electrons, then the quasiparticle 
band opens a gap at the Fermi level. This lowers the 
kinetic energy of all the occupied d-states. The quasi- 
particle band has to be sufficiently narrow for this ferro- 
magnetic charge-density wave to have lower energy than 
the paramagnetically filled quasiparticle band (i.e. the 
gap must be on the order of the quasiparticle bandwidth, 
« t*Z <C t*). Both the charge density fluctuation effect 
and the associated Kondo exchange rapidly vanish away 
from quarter-filling since the charge density wave can no 
longer be commensurate and open a gap at the Fermi 
level. Thus, the ferromagnetism for smaller values of rid is 
due to RKKY exchange (which becomes ferromagnetic). 
Note that this Kondo-effect generated ferromagnetism at 



quarter-filling is a special feature of the bipartite lattice. 

Finally, we speculate that antiferromagnetism should 
return in the region near the empty d-band, since in this 
limit the half-filled f-band PAM reduces to a half-filled 
Hubbard model p"l| ] with strong antiferromagnetic su- 
perexchange [ fL2| . 

Conclusion We investigate the phase diagram and 
screening of the asymmetric PAM in infinite dimensions. 
We find antiferromagnetic, ferromagnetic, and paramag- 
netic ground states, depending upon and other model 
parameters. When the d-band is quarter-filled we find 
a ferromagnetic transition driven by protracted Kondo 
screening when the super and RKKY exchanges are an- 
tiferromagnetic. In the paramagnetic state, we find that 
the temperature dependence of the Kondo screening re- 
tains the same qualitative features as the SIAM, whereas 
the quantitative features are quite different. Specifi- 
cally, when the number of conduction electrons is sig- 
nificantly less than the number of f-electrons, we find 
that Tqp am "C Tqsiam- This results in a protracted re- 
gion of spin screening extending from T > Tqsiam down 
to T < Tqp am or T — T c . In this region, the tem- 
perature dependence of the screened local moment and 
presumably the associated features such as the screen- 
ing resonance are much slower than those of the SIAM, 
which may provide insight into the temperature depen- 
dence of photoemission spectra in some Ce-based heavy- 
fermion lattice compounds j^] . Finally, it is important to 
note that both the suppression of To and the concomitant 
protracted screening disappear when the f-orbital degen- 
eracy N — > oo Jl(|. Thus, Top am also depends strongly 
upon N, so the admixture of crystal-field split states as 
well as a larger N could partially restore an impurity-like 
temperature dependence. 
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FIG. 2. The f-band magnetic moments for the PAM and 
the SIAM vs. temperature for three different values of d-band 
filling when U = 1.5 and V = 0.6. 
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FIG. 3. Half-width of the derivative of the total occupa- 
tion-number, —^f^-, at low temperatures when U = 1.5 and 
V — 0.6. In a Fermi-liquid this quantity should be propor- 
tional to ZT. Any residual value represents the scattering rate 
at the Fermi surface. The flat slope of w(T) at low temper- 
atures for the low-filling case suggests that the Fermi-liquid 
does not begin to form before the ferromagnetic transition (cf. 
Fig. 4). 
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FIG. 1. Kondo temperature vs. d-band filling for the 
infinite-dimensional PAM and SIAM at fixed half-filled f-band 
(rif ?s 1.0) and two different sets of model parameters. The 
inset shows the corresponding Drude weight for the PAM. 
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FIG. 4. Magnetic phase diagram for infinite-dimensional 
PAM at fixed half-filled f-band (n/ « 1.0) and four differ- 
ent sets of model parameters. FM, PM and AFM stand for 
ferromagnetic, paramagnetic and antiferromagnetic phases re- 
spectively. (The FM transition points for U = 2.5, V = 0.5 
are not shown.) The inset shows the difference between the 
staggered and uniform RKKY exchange constants vs. d-band 
filling. 
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